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INTRODUCTION 

This  report  is  the  second  interim  report  outlining  accomplishments 
on  a  research  project  aimed  at  reduction  of  the  earthquake- induced 
liquefaction  hazard  at  Naval  Shore  Facilities.  The  previous  report  (Ref 
1)  provided  a  basic  background  on  the  liquefaction  problem  with  respect 
to  Navy  requirements  and  summarized  the  goals  and  directions  of  research 
to  be  pursued.  As  stated  in  Reference  1,  this  study  was  directed  largely 
toward  monitoring  the  overall  field  of  liquefaction  analysis  in  order 
that  the  most  appropriate  emerging  techniques  be  recognized  for  further 
development  or  adaption  to  Navy  purposes. 

Two  main  aspects  of  liquefaction  development  are  being  considered: 
(1)  site  evaluation  technology  and  (2)  analytical  techniques. 

Reference  2  was  prepared  and  distributed  as  part  of  this  project  to 
serve  as  an  interim  guide  summarizing  current  engineering  practice  for 
assessing  liquefaction  potential.  Major  deficiencies  exist  in  current 
technology,  particularly  in  determining  the  occurrence  of  liquefaction 
in  other  than  level  ground  (no  structure)  situations.  Further,  the 
procedures  for  estimation  of  structural  behavior  -  given  the  occurrence 
of  site  liquefaction  -  are  not  well  established.  Thus,  the  focus  of 
this  effort  is  twofold:  (1)  to  evaluate  the  prediction  of  the  occurrence 
of  liquefaction  by  the  best  field  and  laboratory  techniques,  and  (2) 
given  the  occurrence  of  an  increase  in  pore  pressure,  to  evaluate  its 
effect  on  the  structure.  To  accomplish  this  latter  goal,  emphasis  has 
been  placed  on  developing  effective  stress  analytical  material  models. 

Field  work  has  continued  at  the  Naval  Air  Station  (NAS)  North 
Island,  including  electric  friction  cone  probings  and  piezometric  cone 
soundings;  and  in-house  cyclic  triaxial  testing  of  North  Island  soils 
has  commenced.  These  items  are  further  discussed  in  later  sections  of 
this  report. 


TECHNOLOGY  SUMMARY 

During  the  past  several  years  new  terminology  has  developed  with 
respect  to  liquefaction  of  soils,  and  this  has  led  to  the  potential  for 
misunderstanding  and  misinterpretation.  In  order  to  clarify  this  situa¬ 
tion,  the  Committee  on  Soil  Dynamics  of  the  Geotechnical  Division  of  the 
American  Society  of  Civil  Engineers  has  prepared  a  list  of  recommended 
definitions  (Ref  3). 

The  suggested  terminology  applies  the  term  "liquefaction"  to  the 
situation  wherein  a  soil  is  transformed  into  a  transient  liquefied 
state,  regardless  of  the  initiating  disturbance.  The  use  of  modifiers 
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such  as  "cyclic"  liquefaction,  "initial"  liquefaction,  etc.,  are  dis¬ 
couraged.  The  situation  of  progressively  increasing  strain- to-stress 
ratios  that  occurs  during  cyclic  loading  is  referred  to  as  "cyclic 
strain  softening." 

A  generated  pore  pressure  increment  equal  to  the  effective  con¬ 
fining  pressure  that  existed  prior  to  commencement  of  dynamic  loading  is 
referred  to  as  "full  or  100%  pore  pressure  ratio"  (rather  than  "initial 
liquefaction,"  as  was  used  formerly). 

Distinction  is  made  between  "shear  strength”  and  "shear  resis¬ 
tance."  For  example,  shear  strength  is  defined  as  the  maximum  resistance 
of  soil  to  shearing  stresses.  "Cyclic  shear  resistance"  is  the  level  of 
cyclic  stress  required  to  produce  a  given  strain  (or  liquefaction)  in  a 
specified  number  of  load  cycles. 

"Limited  flow  strain"  or  "limited  flow  deformation"  occurs  follow¬ 
ing  transient  liquefaction,  p/*i^>r  to  recovery  of  resistance  due  to 
dilatency.  "Unlimited  flow  strain"  or  "unlimited  flow  deformations" 
continue  unabated  under  undrained  stress  and  are  accompanied  by  perma¬ 
nent  loss  of  shear  strength. 

Traditionally,  liquefaction  has  been  defined  as  the  situation 
involving  unrecoverable  loss  of  resistance  to  deformation  caused  by 
buildup  of  high  pore  water  pressures. 

According  to  Castro  (Refs  A, 5),  any  soil  whose  ultimate  undrained 
residual  strength  is  less  than  the  in  situ  shear  stresses  is  unstable. 

The  undrained  shear  strength  is  assumed  to  be  a  function  of  only  the 
initial  soil  state  and  independent  of  the  way  in  which  the  failure  load 
is  applied.  This  means  that  the  unlimited  flow  threat  can  be  evaluated 
in  terms  of  two  factors:  (1)  the  undrained  residual  shear  strength  and 
(2)  the  in  situ  stress  state. 

Thus,  the  major  difficulty  in  analyzing  cohesionless  soils  with 
respect  to  earthquake  loading  would  appear  to  be  in  the  prediction  of 
strain  accumulation;  i.e.,  the  cyclic  shear  resistance. 

As  has  been  noted  by  Finn  (Ref  6),  unless  stress  reversals  occur  in 
a  cyclically  loaded  soil  specimen,  pore  pressure  buildup  may  never  rise 
to  the  level  of  confining  pressure  (formerly  called  initial  liquefaction). 
Nevertheless,  very  large  accumulation  of  strain  may  occur  under  stress 
levels  that  are  below  the  static  shear  strength. 

In  light  of  the  above,  experimental  evaluations  of  the  liquefaction 
hazard  appear  to  fall  into  two  areas.  The  first  area  is  theoretically 
trivial,  in  that  any  situation  wherein  field  stress  exceeds  residual 
undrained  strength  is  unstable  and  cannot  be  tolerated.  (This  still 
offers  problems  in  practice  in  that  a  valid  evaluation  of  the  residual 
strength  and  in  situ  shear  stress  may  be  difficult.)  The  second  area, 
that  of  assessing  strain,  becomes  almost  intractable  in  that  as  soon  as 
large  strain  accumulations  commence  in  a  laboratory  specimen,  the  validity 
of  the  test  results  becomes  subject  to  question  (Refs  4,7). 

The  foregoing  is  not  offered  as  an  argument  against  cyclic  testing 
of  soils  in  the  laboratory,  but  rather  to  emphasize  the  degree  of  judgment 
that  is  still  necessary  in  geotechnical  earthquake  engineering.  This 
judgment  becomes  even  more  important  with  regard  to  complex  waterfront 
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construction.  For  example,  it  is  shown  (Ref  8)  that  earth  dams,  which 
are  relatively  forgiving  as  far  as  deformation  is  concerned,  generally 
stand  up  to  earthquakes  quite  well.  Thus,  only  in  cases  of  poorly 
constructed,  saturated,  cohesionless  materials  or  under  unusually  severe 
ground  motions  do  earth  dams  generally  undergo  major  earthquake  damage. 
However,  soil  retaining  structures,  such  as  are  prevalent  at  the  water¬ 
front,  have  wide  ranges  of  deformation  acceptability.  Some  structures, 
such  as  sheet  pile  cells  or  drydock  walls,  stay  fail  under  relatively 
minor  dislocations,  whereas  relieving  platforms,  tied  bulkhead  walls, 
etc.,  in  some  cases  may  permit  relatively  large  deformations  prior  to 
functional  failure. 

Another  complexity  is  that  the  most  dramatic  earthquake-related 
liquefaction  failure  may  occur  well  after  the  earthquake  motions  have 
ended,  such  as  in  Ojiha  (1939)  or  Niigata  (1964),  Japan  (Ref  4). 

The  influence  of  the  degree  of  nonsymmetry  of  load  application  does 
not  appear  to  have  been  addressed  in  any  general  manner,  but  rather  the 
test  data  have  been  applied  directly  to  specific  cases. 

Past  engineering  practice  generally  has  been  to  design  geotechnical 
structures  to  withstand  earthquakes  by  including  the  equivalent  static 
inertial  forces  in  the  traditional  stability  analysis  and  checking  to 
see  that  the  factor  of  safety  is  sufficiently  above  unity.  Alternative 
static  procedures  using  slip  planes  in  the  soil  mass  compute  acceleration 
levels  that  may  cause  yield.  Seed  (Ref  8)  concludes  that  pseudo-static 
analysis  techniques  must  be  used  with  great  caution  and  that  dynamic 
analysis  techniques  provide  a  more  reliable  basis  for  estimating  perfor¬ 
mance  and  safety.  The  critical  period  of  stability  may  occur  due  to 
redistribution  of  pore  water  pressures  after  the  earthquake  shaking  has 
ceased.  This  type  of  failure  cannot  be  addressed  by  pseudo-static 
analysis. 

Evaluation  of  saturated  cohesionless  soils  can  only  properly  be 
made  when  pore  pressures  and  resulting  reduced  confining  stress  and 
stiffness  changes  are  taken  into  account.  Cohesive  soils  also  exhibit 
complex  behavior.  Thus,  large  deformations  can  occur  under  prolonged 
oscillating  loads  even  though  maximum  applied  stresses  are  less  than  the 
static  strength  of  the  soil. 

Analysis  of  total  stress  does  not  directly  yield  information  on 
pore  pressures  and,  hence,  strain  softening,  as  the  reduction  in  effective 
stress  reduces  shear  modulus.  Only  a  nonlinear  model  based  upon  effective 
stress  can  fully  couple  pore  pressure  generation  and  softening,  which 
can  significantly  affect  the  dynamic  response.  Drainage  of  pore  pressure 
postulated  into  effective  stress  models  will  be  required  to  complete  the 
representation  of  a  soil. 

Significant  work  certain  to  have  a  major  impact  on  analysis  tech¬ 
niques  is  in  progress  in  the  development  of  effective  stress  material 
models  for  use  with  finite  element  codes.  The  effective  stress  models 
will  calculate  actual  dynamically  induced  pore  pressures.  Work  is  in 
progress  by  several  researchers:  Finn,  Leek,  and  Martin  (Ref  9); 

Ghaboussi  and  Dikman  (Ref  10);  Baladi  and  Rohani  (Ref  11);  Mroz,  Norris, 
and  Zienkiewicz  (Ref  12);  and  Prevost  and  Hughes  (Ref  13). 


The  nonlinear  effective  stress  material  models  will  require  detailed 
material  response  data  to  fit  the  parameters  required.  Once  an  effective 
stress  material  model  has  been  implemented  in  a  general  finite  element 
code  and  verified  by  the  engineering  community,  it  could  become  the  most 
important  tool  in  liquefaction  analysis. 


FIELD  INVESTIGATIONS 

During  the  investigation  of  the  earthquake-induced  liquefaction 
threat  at  NAS  North  Island  (Ref  14),  several  unusual  soils  were  encoun¬ 
tered.  Portions  of  NAS  originally  below  sea  level  had  been  reclaimed, 
using  hydraulic  fills  -  primarily  silty  sands  -  dredged  from  the  bottom 
of  San  Diego  Bay.  These  reclaimed  areas  had  been  subjected  to  several 
series  of  fillings,  most  between  1919  and  1952.  As  a  result  of  investi¬ 
gations  made  on  North  Island  in  connection  with  foundations  of  existing 
structures,  a  recently  filled  region  was  selected  for  further  study  (Ref 
1).  Available  soil  data  were  supplemented  by  additional  soil  investiga¬ 
tions,  including  dynamic  split-spoon  penetration  tests  and  static  cone 
penetrometer  soundings.  The  split-spoon  penetrations  were  conducted  as 
specified  in  ASTM  D  1586  except  that  a  2-1/2-inch  (California-type) 
split  spoon  was  used.  The  cone  penetration  soundings  conformed  to  ASTM 
D  3441-75T.  Based  upon  this  information,  the  generalized  soil  profile 
shown  in  Figure  1  was  developed.  This  soil  profile  encompasses  an  old 
filled  bay  channel,  formerly  known  as  Spanish  Bight,  which  was  hydrauli¬ 
cally  filled  to  its  present  elevation  in  1945.  The  old  bay  bottom 
immediately  prior  to  the  first  filling  is  represented  by  layer  5. 

The  results  of  split-spoon  penetration  tests  and  friction  cone 
soundings,  denoted  as  penetration  holes  PI  through  P5  along  the  soil 
profile  in  Figure  1,  are  shown  in  Figure  2.  The  split-spoon  penetration 
values  near  the  boundary  between  layers  4  and  5  fall  as  low  as  one  or 
two  blows  per  foot,  suggesting  a  relative  density,  D  *,  for  soil  stratum 
5  of  less  than  50%  (see  Ref  14).  r 

The  relatively  low  friction  cone  penetration  readings  of  5  to  10 
kg/sq  cm  shown  in  Figure  2  suggest  values  of  relative  density  on  the 
order  of  only  20%  or  30%.  As  suggested  in  Reference  14,  natural  soils 
with  relative  densities  of  less  than  40%  are  expected  to  be  quite  rare 
in  nature.  However,  based  upon  the  available  information  it  was  necessary 
to  assume,  for  analysis,  that  the  soil  in  the  region  of  layer  5  could 
not  have  a  relative  density  greater  than  35%. 

This  type  of  soil  appears  in  waterfront  areas  throughout  California 
(such  as  at  Long  Beach  Naval  Shipyard)  and,  thus,  is  of  special  interest 
to  the  Navy.  Most  regions  in  California  can  be  expected  to  experience 
ground  motions  of  some  severity  relatively  often  (in  the  geologic  time 
scale).  Therefore,  it  was  considered  critical  that  the  liquefaction 
potential  of  this  soil  be  explored  further. 


*Relative  density  is  a  measure  of  the  present  density  of  a  soil  with 
respect  to  its  maximum  and  minimum  densities. 
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A  program  of  undisturbed  sampling  and  laboratory  cyclic  triaxial 
testing  was  carried  out.  The  sampling  program  consisted  of  obtaining 
soil  samples  along  the  soil  profile  shown  in  Figure  1,  using  an  Osterberg 
piston  sampler.  A  description  of  this  sampling  and  testing  program  is 
presented  in  Reference  14  and  discussed  elsewhere  in  this  report. 

This  section  of  the  report  addresses  improved  reconnaissance  or 
field  liquefaction  threat  evaluation  techniques.  No  matter  how  sophis¬ 
ticated  and  reliable  analytical  procedures  may  become,  they  can  provide 
accurate  response  predictions  only  when  valid  input  data  can  be  obtained. 

In  the  case  of  minor  structures,  where  involved  analytical  programs  are 
not  possible,  expedient  field  evaluation  techniques  become  the  only 
available  approach. 

In  this  regard,  friction  cone  soundings  using  a  different  type  of 
cone  were  performed  during  April  1979  at  North  Island  in  the  vicinity  of 
P3  on  the  soil  profile  of  Figure  1  (see  Figure  3).  This  investigation 
used  an  electric  friction  cone  rather  than  a  mechanical  cone  as  was  used 
to  provide  the  data  in  Figure  2.  The  electric  cone  has  a  standard 
10-sq-cm,  60-degree  tip  and  a  150-sq-cm  sleeve,  like  the  mechanical 
cone.  However,  rather  than  measuring  the  point  and  sleeve  resistances 
over  10-em  depth  increments,  as  with  the  mechanical  cone,  the  electric 
cone  can  record  penetration  resistances  continuously .  Penetration  must 
be  interrupted  only  at  1-meter  increments  to  add  additional  rods.  A 
schematic  picture  of  the  electric  friction  cone  is  shown  in  Figure  4.  A 
total  of  nine  friction  cone  soundings  was  conducted.  Complete  results 
have  been  presented  in  Reference  15;  however,  Figure  5  shows  a  typical 
penetration  resistance  record.  The  electric  cone  records  were  found  to 
be  very  similar  to  the  mechanical  cone  penetration  records  illustrated 
for  this  area  (P3  in  Figure  1).  Both  records  showed  very  stiff  layers 
near  depths  of  3  and  16  feet  and  a  very  soft  clayey  (high-friction 
ratio)  layer  near  a  depth  of  13  feet. 

The  liquefaction  potential  of  a  soil  is  directly  related  to  its 
volumetric  change  tendencies,  more  specifically  to  its  volumetric-strain/ 
shear-strain  coupling.  Traditionally,  the  volume  change  properties  have 
been  the  least  studied  from  the  standpoint  of  in  situ  measurements. 

Problems  associated  with  sample  disturbance,  laboratory  simulation  of  in 
situ  stress  state,  temperature,  chemical  and  biological  environments, 
and  soil  heterogeneity,  for  example,  severely  impair  any  analysis  based 
upon  laboratory  testing. 

Thus,  although  in  situ  testing  is  attractive  in  theory,  direct 
determination  of  in  situ  volume  change  properties  applicable  to  earthquake- 
type  loading  is  not  yet  possible.  Recourse  must  be  made  to  empirical 
correlations  between  volume  change  properties  and  some  form  of  expedient 
field  test.  Various  forms  of  penetration  tests  have  been  utilized  in 
this  regard,  particularly  the  standard  penetration  test  (ASTM  D  1586) 
and  the  friction  cone  penetration  test  (ASTM  D  3441-75T). 

Unfortunately,  the  penetration  resistance  of  a  soil  is  a  function 
not  only  of  its  volume  change  characteristics  but  also  of  its  strength, 
shear  stiffness,  and  other  deformational  characteristics.  Penetration 
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resistance  may  be  influenced  by  factors  other  than  those  directly  influ¬ 
encing  liquefaction  potential.  Therefore,  it  may  be  necessary  to  measure 
several  types  of  response  before  accurate  penetration  test  correlations 
become  possible. 

One  promising  device  with  regard  to  in  situ  volume  change  is  a 
piezometer  probe  (Ref  16)  which  measures  pore  pressure  in  the  tip  of  a 
penetrating  cone.  If,  during  cone  penetration,  positive  (increased) 
pore  water  pressures  are  generated,  then  the  effective  stresses  and, 
hence,  the  strength  and  resistance  to  penetration  are  reduced.  Alter¬ 
natively,  if  negative  (reduced)  pore  water  pressures  occur,  then  the 
soil  structure  is  dilating,  and  effective  stresses  are  increased  and  so 
is  penetration  resistance.  The  measured  incremental  changes  in  pore 
pressure  during  penetration  of  saturated  soils  are  directly  related  to 
volume  change  tendency.  They  are  also  a  function  of  soil  permeability 
and  rate  of  penetration. 

The  rate  at  which  pore  water  pressures  reach  equilibrium  following 
cessation  of  penetration  are  direct  functions  of  permeability.  Thus,  a 
device  exists  with  the  potential  for  correlation  with  both  volume  change 
characteristics  and  permeability.  These  two  factors  (along  with  nature 
of  loading)  are  the  major  determinants  which  control  the  occurrence  and 
severity  of  soil  liquefaction. 

In  order  to  better  appreciate  the  potential  of  the  piezometric  cone 
or  probe  for  liquefaction  prediction,  a  program  of  piezometric  cone 
penetrations  was  inititated,  through  contract,  at  NAS  North  Island  (Ref 
15).  Eight  soundings  were  performed  in  the  vicinity  of  holes  P2  and  P3 
(see  Figures  1  and  3).  This  work  was  accomplished  by  Fugro  Gulf,  Inc. 
of  Houston,  Tex.,  using  a  truck-mounted  electronic  cone  penetrometer 
system.  The  piezometric  cone  is  shown  in  Figure  6.  The  cone  tip  has 
the  same  dimensions  as  the  standard  friction  cone.  However,  a  specially 
designed  porous  element  is  used  to  allow  entry  of  pore  water.  Thus, 
both  pore  pressure  and  total  tip  resistance  are  measured.  To  simplify 
the  instrumentation,  sleeve  friction  is  not  measured  with  this  pene¬ 
trometer.  Pore  pressure  is  recorded  as  a  function  of  depth  and  penetra¬ 
tion  rate  during  penetration,  and  pore  pressure  dissipation  is  recorded 
as  a  function  of  time  when  penetration  is  stopped. 

In  operation,  special  precautions  were  required  to  prevent  air 
entry  (water  cavitation)  into  the  porous  element  during  cone  passage 
through  the  5-foot-thick  soil  zone  above  the  water  table.  These  pre¬ 
cautions  included  pre-drilling  and  casing  holes  to  below  the  water 
table.  The  piezometric  cone  element  was  then  saturated  and  encased  in  a 
water-filled  membrane  which  maintained  tip  saturation  while  the  pene¬ 
trometer  was  being  lowered  in  the  bored  hole  to  beneath  the  water  table. 
The  membrane  was  subsequently  punctured  by  the  advancing  tip  when  pene¬ 
tration  commenced  into  the  soil  at  the  bottom  of  the  pre-drilled  hole. 

Records  of  the  eight  piezometric  cone  soundings  are  too  extensive 
to  be  presented  herein,  but  a  typical  data  plot  (probe  4A)  is  shown  in 
Figure  7.  A  complete  record  is  presented  in  Reference  15. 

The  pore  pressures  during  penetration  indicate  both  dilative 
(decreased  pore  pressure)  behavior  and  compactive  (pore  pressure  increase) 
behavior  for  the  different  in  situ  soil  strata. 
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Figure  7  indicates  a  medium-dense  layer  of  soil  between  about  6  and 
10  feet  < eep  which  tends  to  dilate  due  to  passage  of  the  cone  and  gener¬ 
ates  a  reduction  in  pore  pressures.  Beyond  a  depth  of  10  feet,  the  soil 
strength  falls  off,  with  attendant  pore  pressure  generation.  This 
corresponds  to  the  very  sensitive  silts  (layer  5)  in  soil  profile  of 
Figure  1.  A  denser  stratum  is  again  encountered  at  a  depth  of  about  15 
feet  (approaching  layer  6)  with  a  dramatic  reduction  in  pore  pressure 
generation.  Tip  resistance  then  falls  off  due  to  no  further  penetration. 

Thus,  it  appears  that  the  piezometric  cone  is  sensitive  to  the  type 
of  soil  response  associated  with  liquefaction.  Whether  the  results  of 
such  studies  are  precise  enough  to  be  used  for  reliable  liquefaction 
potential  assessments  cannot  be  ascertained  until  further  experience  has 
been  acquired. 

LABORATORY  LIQUEFACTION  RESISTANCE 

The  most  widely  utilized  laboratory  technique  for  evaluating  the 
liquefaction  resistance  of  soils  is  based  upon  stress-controlled,  cyclic 
triaxial  testing.  Liquefaction  potential  was  measured  for  various  can¬ 
didate  soil  types  using  cyclic  triaxial  equipment  at  CEL.  A  description 
of  this  equipment  together  with  representative  test  results  for  a  "stan¬ 
dard"  soil  has  been  presented  in  Reference  17.  Reference  17  would  also 
serve  as  a  guide  for  performing  cyclic  triaxial  testing  to  be  used  in 
conjunction  with  evaluating  the  seismic  liquefaction  potential  at  Naval 
facilities.  This  test  equipment  has  been  used  for  other  CEL  projects 
dealing  with  runway  subgrade  evaluations  and  resilient  modular  testing 
of  airfield  base  course  materials.  Its  major  application,  however,  is 
for  earthquake-related  soil  liquefaction  studies,  such  as  those  conducted 
at  NAS  North  Island  (see  Ref  1  and  14). 

An  in-house  testing  program  has  been  initiated  to  investigate  the 
liquefaction  response  of  typical  soils  encountered  at  the  waterfront 
and,  thus,  of  special  interest  to  the  Navy.  During  the  investigation 
outlined  in  Reference  14,  a  deposit  of  sensitive  silty  sand  (layer  5  in 
Figure  1)  was  encountered  which  had  extremely  low  penetration  resistance. 

A  limited  number  of  cyclic  triaxial  tests  on  this  soil  were  performed 
under  contract,  as  outlined  in  Reference  1.  Additional  specimens  acquired 
during  this  investigation  were  selected  for  further  in-house  cyclic 
triaxial  testing. 

The  soil  selected  for  cyclic  testing  was  a  very  fine  silty  sand, 
generally  dark  grey  in  color.  A  typical  grain  size  analysis  is  shown  in 
Figure  8.  This  soil  is  relatively  uniform  with  a  coefficient  of  uni¬ 
formity  C  of  less  than  4.0.  Maximum  and  minimum  dry  densities,  deter¬ 
mined  as  outlined  in  Reference  18  were  108.2  pcf  and  79.4  pcf,  respec¬ 
tively.  The  density  of  undisturbed  samples  was  approximately  93  pcf, 
which  was  equivalent  to  a  relative  density  of  about  55%. 

Two  series  of  cyclic  triaxial  tests  were  performed.  Series  1 
consisted  of  seven  tests  performed  on  undisturbed  tube  samples  2.88 
inches  in  diameter  with  a  length-to-diameter  ratio  of  approximately  2. 
Average  dry  densities  varied  from  90.8  to  94.2  pcf  with  an  average 
density  of  92.5  pcf  (D^  =  53%). 


-  -  ■■■>■ .  - .  . 


Series  "  consisted  of  five  cyclic  tests  on  specimens  2.80  inches  in 
diameter  and  6.0  inches  long.  These  specimens  were  reconstituted  from 
the  Series  1  samples  by  means  of  wet  tamping. 

All  test  specimens  were  consolidated  isotropically  under  effective 
confining  pressures  of  6  to  8  psi,  corresponding  to  the  former  in  situ 
vertical  effective  stress  value.  Harmonic  cyclic  shear  loading  was 
applied  at  a  frequency  of  1  hertz,  and  the  number  of  cycles  required  to 
reach  full  pore  pressure  ratio  and  5%  double  amplitude  strain  was  recorded. 
Pulsating  deviator  stress,  axial  strain,  and  pore  water  pressure  were 
recorded  on  a  Honeywell  Visicorder  strip-chart  recorder. 

The  test  results  for  Series  1  are  shown  in  Figures  9  and  10.  Also 
shown  on  Figures  9  and  10  are  the  best  fit  curves  for  the  results  of  the 
cyclic  tests  reported  in  Reference  1  (conducted  by  a  contractor).  Test 
results  for  the  Series  1  and  2  tests  are  compared  in  Figures  11  and  12. 

As  indicated  in  Figures  9  and  10,  the  cyclic  strength  curves  from 
CEL  show  strengths  between  30%  and  45%  stronger  than  those  obtained  by 
the  contractor.  Although  this  might  be  explained  by  the  limited  number 
of  samples  tested,  the  relative  consistency  or  repeatability  of  the  two 
separate  data  sources  tend  to  preclude  this.  The  possibility  was  con¬ 
sidered  that  the  longer  period  of  storage  of  the  CEL-tested  samples 
resulted  in  strengths  higher  than  the  contractor's  test  results.  This 
latter  explanation  is  difficult  to  accept  in  view  of  the  very  minor 
differences  shown  between  the  CEL-tested  undisturbed  and  reconstituted 
specimen  test  results  shown  in  Figures  11  and  12.  However,  the  results 
shown  in  Figures  11  and  12  might  be  influenced  by  specimen  inhomogenei¬ 
ties  . 

For  example,  the  reconstituted  specimens  formed  by  moist  tamping 
would  be  expected  to  be  relatively  homogeneous  and,  therefore,  provide  a 
relatively  uniform  strength  throughout.  The  undisturbed  sample  strengths 
could  be  influenced  by  small  weaker  zones  which  are  the  result  of  envi¬ 
ronmental  variations  during  the  initial  sedimentation  process.  Such 
weaker  zones  would  tend  to  counter  the  tendency  of  undisturbed  samples 
to  have  higher  strengths  than  reconstituted  samples.  Nevertheless,  even 
though  the  differences  between  the  in-house  and  the  contractor  test 
results  could  conceivably  be  ascribed  to  time-dependent  structural 
differences,  Reference  18  concludes  that  these  differences  are  largely 
due  to  variations  in  test  procedure  and  equipment.  Thus,  it  appears 
that  with  these  sensitive  materials,  carefully  conducted  experimental 
testing  programs  at  two  different  facilities  resulted  in  significant 
differences  in  liquefaction  related  strength  measurements.  A  standar¬ 
dized  procedure  is  essential  for  consistent  results.  References  17  and 
18,  prepared  as  part  of  this  task,  assist  in  this  area. 


EFFECTIVE  STRESS  SOIL  MODELS 

Current  engineering  practice  does  not  allow  for  the  accurate  esti¬ 
mating  of  displacements  of  a  structure  founded  on  subsurface  layers  of 
potentially  liquefiable  soil.  To  accomplish  such  estimates,  effective 
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Series  2  consisted  of  five  cyclic  tests  on  specimens  2.80  inches  in 
diameter  and  6.0  inches  long.  These  specimens  were  reconstituted  from 
the  Series  1  samples  by  means  of  wet  tamping. 

All  test  specimens  were  consolidated  isotropically  under  effective 
confining  pressures  of  6  to  8  psi,  corresponding  to  the  former  in  situ 
vertical  effective  stress  value.  Harmonic  cyclic  shear  loading  was 
applied  at  a  frequency  of  1  hertz,  and  the  number  of  cycles  required  to 
reach  full  pore  pressure  ratio  and  5%  double  amplitude  strain  was  recorded. 
Pulsating  deviator  stress,  axial  strain,  and  pore  water  pressure  were 
recorded  on  a  Honeywell  Visicorder  strip-chart  recorder. 

The  test  results  for  Series  1  are  shown  in  Figures  9  and  10.  Also 
shown  on  Figures  9  and  10  are  the  best  fit  curves  for  the  results  of  the 
cyclic  tests  reported  in  Reference  1  (conducted  by  a  contractor).  Test 
results  for  the  Series  1  and  2  tests  are  compared  in  Figures  11  and  12. 

As  indicated  in  Figures  9  and  10,  the  cyclic  strength  curves  from 
CEL  show  strengths  between  30%  and  45%  stronger  than  those  obtained  by 
the  contractor.  Although  this  might  be  explained  by  the  limited  number 
of  samples  tested,  the  relative  consistency  or  repeatability  of  the  two 
separate  data  sources  tend  to  preclude  this.  The  possibility  was  con¬ 
sidered  that  the  longer  period  of  storage  of  the  CEL-tested  samples 
resulted  in  strengths  higher  than  the  contractor's  test  results.  This 
latter  explanation  is  difficult  to  accept  in  view  of  the  very  minor 
differences  shown  between  the  CEL-tested  undisturbed  and  reconstituted 
specimen  test  results  shown  in  Figures  11  and  12.  However,  the  results 
shown  in  Figures  11  and  12  might  be  influenced  by  specimen  inhomogenei¬ 
ties  . 

For  example,  the  reconstituted  specimens  formed  by  moist  tamping 
would  be  expected  to  be  relatively  homogeneous  and,  therefore,  provide  a 
relatively  uniform  strength  throughout.  The  undisturbed  sample  strengths 
could  be  influenced  by  small  weaker  zones  which  are  the  result  of  envi¬ 
ronmental  variations  during  the  initial  sedimentation  process.  Such 
weaker  zones  would  tend  to  counter  the  tendency  of  undisturbed  samples 
to  have  higher  strengths  than  reconstituted  samples.  Nevertheless,  even 
though  the  differences  between  the  in-house  and  the  contractor  test 
results  could  conceivably  be  ascribed  to  time-dependent  structural 
differences,  Reference  18  concludes  that  these  differences  are  largely 
due  to  variations  in  test  procedure  and  equipment.  Thus,  it  appears 
that  with  these  sensitive  materials,  carefully  conducted  experimental 
testing  programs  at  two  different  facilities  resulted  in  significant 
differences  in  liquefaction  related  strength  measurements.  A  standar¬ 
dized  procedure  is  essential  for  consistent  results.  References  17  and 
18,  prepared  as  part  of  this  task,  assist  in  this  area. 


EFFECTIVE  STRESS  SOIL  MODELS 

Current  engineering  practice  does  not  allow  for  the  accurate  esti 
mating  of  displacements  of  a  structure  founded  on  subsurface  layers  of 
potentially  liquefiable  soil.  To  accomplish  such  estimates,  effective 
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stress  soil  models  are  being  evaluated.  Two  paths  are  being  investigated 
In  the  first,  a  full  nonlinear  elasto-plastic  model  of  soil  skeleton 
behavior  is  desired  which  accounts  for  nonlinear  behavior  and  volume 
changes  with  load.  A  second  path  takes  a  more  simplified  approach, 
separating  the  basic  material  elasto-plastic  behavior  and  the  cumulative 
volume  change.  This  section  will  summarize  results  on  the  two  approaches 
implemented  by  three  soil  models. 

Baladi  Plasticity  Model 

Study  of  an  elasto-plastic  cap  model  for  cohesionless  soil  is  being 
conducted  at  the  Army  Waterways  Experiment  Station.  A  minor  part  of  the 
funding  was  provided  by  this  work  unit.  A  detailed  description  of  the 
model  and  supporting  equations  are  provided  in  Reference  19.  The  model 
is  characterized  by  seven  parameters  which  define  bulk  modulus  variation 
with  the  first  invariant  of  the  stress  tensor,  shear  modulus  variation 
with  second  invariant  plastic  strain,  the  failure  envelope,  and  shape  of 
the  cap.  Additional  parameters  can  be  introduced  for  increased  defini¬ 
tion  of  specific  soil  behavior.  Material  properties  for  the  model  are 
derived  from  drained  hydrostatic  and  triaxial  tests.  Present  development 
is  only  partially  complete  since  the  model  computes  effective  stresses, 
but  pore  pressures  are  not  provided  for.  Presently,  the  model  is  set  up 
for  explicit  calculation,  and  implicit  procedures  would  also  be  desirable 
Figure  13  shows  typical  computed  undrained  results,  showing  stress 
difference  as  a  function  of  axial  strain.  Figure  13a  and  b  shows  the 
stress  difference  versus  axial  strain  and  stress  pressure,  respectively 
for  a  biaxial  test  on  sand.  When  the  soil  reaches  the  failure  envelope, 
it  begins  to  fail  and  the  loading  stress  difference  falls  as  noted  in 
Figure  13b.  Figure  13c  shows  the  stress  difference-axial  strain  for 
another  soil.  Figure  13d  shows  its  stress  path.  In  this  case  the  soil 
increases  in  strength.  The  characterization  of  the  soil  behavior  is 
controlled  by  the  material  parameters  input  to  the  model. 

Zienkiewicz  Model 


This  model  is  a  more  pragmatic  approach.  The  model  consists  of  a 
simple  ideally  elasto-plastic  model  for  the  soil  with  a  nonassociative 
flow  rule  corresponding  to  the  Mohr-Coulomb  surface.  A  separate  ex¬ 
pression  for  determining  progressive  increase  in  volume  strain  is  added. 
The  formulation  is  valid  for  undrained  conditions.  Data  based  on  simple 
shear  tests  are  used  as  the  basis  for  developing  a  relationship  between 
volume  strain  and  a  strain-related  damage  parameter.  Additional  detail 
may  be  found  in  Reference  20.  This  model  is  presently  being  implemented 
and  evaluated. 

Prevost  Model 


An  anistropic  nonlinear  elasto-plastic  path-dependent  model  for 
soils  subject  to  complicated  loading  has  been  developed  by  Prevost. 
Funding  was  provided  by  CEL  for  implementation  of  this  model  as  part  of 
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this  task.  An  elliptic  yield  surface  is  used  in  conjunction  with  critical 
state  lines  and  an  associative  flow  rule  to  compute  plastic  strain  rate 
vectors.  A  nonassociative  flow  rule  is  then  used  on  the  yield  surface. 
Material  parameters  specify  the  size  and  position  of  the  yield  surfaces 
with  associated  moduli.  Table  1  and  Figures  14  and  15  give  results  for 
Cook's  Bayou  sand.  Additional  information  may  be  found  in  Reference  21. 
This  model  is  under  evaluation  at  the  present  time. 


OVERVIEW 

Available  procedures  for  analyzing  the  liquefaction  hazard  to  Navy 
facilities  are  inadequate.  With  conventional  approaches  to  predict 
liquefaction  at  NAS  North  Island,  factors  of  safety  under  earthquake 
loadings  appear  to  vary  by  more  than  200%.  What  is  needed  is  a  redefin¬ 
ing  of  the  liquefaction  phenomenon  and  development  of  less  empirical 
treatments . 

In  this  regard,  the  loss  of  strength  of  a  saturated  cohesionless 
soil  under  dynamic  loading  can  be  treated  under  two  types  of  response: 

(1)  the  inherently  unstable  situation  wherein  in  situ  shear  stresses 
exceed  undrained  residual  strengths,  and  (2)  the  case  of  accumulative 
displacements.  In  the  latter  case,  effective  stresses  either  are  not 
reduced  substantially  during  dynamic  loading,  or  are  reduced  for  such 
short  durations  that  unrestrained  failures  do  not  occur.  Nevertheless, 
deformations  may  accumulate  during  intervals  of  reduced  soil  stiffness 
and  result  in  unacceptable  levels  of  displacement. 

It  is  concluded  that  the  inherently  unstable  situation  of  the  first 
case  is  unacceptable.  Any  situation  wherein  strain  softening  could 
result  in  residual  strengths  below  applied  stresses  requires  that  remedial 
action  be  taken  or,  at  least,  that  allowance  be  made  for  outright  failure. 
The  case  of  accumulation  of  displacements  requires  considerable  additional 
effort.  Unless  large  scale  testing  techniques  capable  of  reliably 
modeling  dynamic  soil  deformations  throughout  extended  ground  motion 
intervals  are  made  available,  the  only  valid  avenue  of  approach  would 
seem  to  be  through  effective  stress  soil  models.  Short  of  instrumentation 
and  testing  of  prototype  structures,  only  analytical  models  promise  the 
capability  of  handling  pore  water  migrations  and  soil  structure  interac¬ 
tion  involved  with  the  liquefaction  problem. 
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Table  1 .  Continued 


SYMBOLS: 

G 

= 

Elastic  shear  modulus 

(coni' d) 

P 

= 

Coordinate  locating  center  of  yield  surface 

o' 

VC 

= 

Effective  confining  stress 

°3c 

= 

Chamber  confining  pressure  in  triaxial  test 

°dp 

= 

Deviator  stress  in  triaxial  test  (dynamic) 

k 

c 

2 

Ratio  of  horizontal  to  vertical  stress  in 
triaxial  test 

E 

V 

= 

Volumetric  strain 

a  ,a 
v’  y 

= 

Stresses  in  coordinate  directions  x,y 

Cone  Resist  a  net- 
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Figure  5.  Cone  penetration  test  no. 
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son  of  cyclic  triaxial  data  at  5 7.  double  amplitude 
between  CEL  and  the  contractor.  (Ref  18). 


Figure  11.  Comparison  of  cyclic  strength  between  undisturbed  field  samples 

and  laboratory  reconstituted  samples  at  initial  liquefaction.  (Ref  18) 


Comparison  of  cyclic  strength  between  undisturbed  field  samples 
and  laboratory  reconstituted  samples  at  5%  double  amplitude  strain 
(Ref  18) . 


Tigure  14.  Calculated  and  mea¬ 
sured  stress-strain  behavior 
for  consolidation/swelling 
tests  with  Cook's  Bavou  Sand. 
(Ref  21). 


(b)  Volumetric  stress-strain  behavior. 


Figure  15.  Calculated  and  mea¬ 
sured  stress-strain  behavior 
for  triaxial  soil  test  with 
Cook's  Bayou  Sand.  (Ref  21) 
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